arXiv:astro-ph/0503384vl 17 Mar 2005 


Astronomy & Astrophysics manuscript no. 2834main February 2, 2008 

(DOI: will be inserted by hand later) 


The first observed stellar X-ray flare oscillation: Constraints on 
the flare loop length and the magnetic field 

U. Mitra-Kraev, L. K. Harra, D. R. Williams, and E. Kraev 

Mullard Space Science Laboratory, University College London, Holmbury St. Mary, Dorking, Surrey RH5 6NT, UK 
Received 7 February 2005 / Accepted 12 March 2005 

Abstract. We present the first X-ray observation of an oscillation during a stellar flare. The flare occurred on the active 
M-type dwarf AT Mic and was observed with XMM-Newton. The soft X-ray light curve (0.2-12 keV) is investigated with 
wavelet analysis. The flare’s extended, flat peak shows clear evidence for a damped oscillation with a period of around 750 s, 
an exponential damping time of around 2000 s, and an initial, relative peak-to-peak amplitude of around 15%. We suggest 
that the oscillation is a standing magneto-acoustic wave tied to the flare loop, and find that the most likely interpretation is a 
longitudinal, slow-mode wave, with a resulting loop length of (2.5 ± 0.2) 10*® cm. The local magnetic field strength is found 
to be 105 ± 50 G. These values are consistent with (oscillation-independent) flare cooling time models and pressure balance 
scaling laws. Such a flare oscillation provides an excellent opportunity to obtain coronal properties like the size of a flare loop 
or the local magnetic field strength for the otherwise spatially-unresolved star. 
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1. Introduction 

Oscillations in the solar corona have by now been observed for 
many years. Wavelength regimes ranging from radio to hard X- 
rays have been investigated in t he sea rch for evidence of waves. 
Table 1 in lAschwanden et’aH ll 19991) provides a summary and 
description of the different periods found (0.02 to 1000 s). Most 
of these waves have been explained by m agneto-hydrody namic 
(MHD) oscillations in coronal loops, see iRobertsI ll200(l) for a 
detailed review of waves and oscillations in the corona. 

One of the most exciting aspects of “coronal seismology” 
is that it potentially provides us with the capabi lity for deter¬ 
minin g the magnetic field strength in the corona jRoberts et alJ 
11984 . as well as, in the stellar case, with information on oth¬ 
erwise unresolved spatial scales, e.g., flare loop lengths. It is 
notoriously difficult to measure the magnetic field strength in 
the corona because of the very high speeds of coronal electrons, 
which broaden spectral lines far beyond the width of Zeeman 
splitting. Techniques using both near-infrared emission lines 
and radio observations are successful, but have poor spatial 
resolution. Indirect methods are commonly used, such as the 
extrapolations of the coronal magnetic field from the photo- 
spheric magnetic field, which in turn can be measured using 
the Zeeman effect. 

Many of the observations of waves have been determined 
from variations in intensity. However, a huge step forward was 
achieved in solar coronal physics due to the high spatial resolu- 
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tion availa ble with the Transit i on Region and Coronal Exp lorer 
(TRACE, iHandv et afl ll 999l) . lAschwanden et all lll999l) ob¬ 
served the first spatial displacement oscillations of coronal 
loops. It was suggested that these oscillations were triggered by 
a fast-mode shock from a flare site, and they were interpreted 
as standing fast kink-mode waves. 

iNakariakov & Ofmanlll200ll) made use of such flare-related 
spatial oscillations to determine the magnetic field strength 
(13 + 9 G), which in the case of a standing kink wave is re¬ 
lated to the period of the oscillation, the density of the loop, 
and the length of the loop. Another model used to derive the 
m agnetic field strength from M op oscillations was put forward 
bv IZaitsev & Stepanovl(ll989l) . which assumes that the oscilla¬ 
tion is triggered by a centrifugal force, generated by the evap¬ 
orating chromospheric plasma moving upward along the mag¬ 
netic field. In this case, the magnetic field strength is given by 
the amplitude of the oscillation, as well as the loop density and 
temperature. The model also predicts loop lengths if collisional 
damping is assumed. (Or vice versa, density and temperature, 
if the loop length is known.) 

Although waves are observed across the electromagnetic 
spectrum on the Sun, observations on other stars are rarer. 
The principal reason for this is that the Sun can be spatially 
resolved, whereas on stars, the signal for the wave must be 
strong enough to be observed above the full disk emission. The 
first oscillation associated with a stellar flare was reported by 
lRodonbllll974 . who observed the flare star II Tau with high¬ 
speed optical photometry, and found a long-lived oscillation 
during flare decay with a mean period of 13 s. In a series of pa- 
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pers in t he late 1980’s a nd early 1990’s, Andrews and cowork¬ 
ers (e.g.. lAndrewsll99'(il) presented optical observations of dMe 
flare stars, where they observed quasi-periodicities with peri¬ 
ods of a few tens of seconds, and put forward the idea that 
these were coronal loop oscillations. iMullan et alJ lll992h also 
found optical oscillations (periods of the order of a few min¬ 
utes) in dMe stars, concluding that they more likely arise from 
a coronal than from a ph otospheric origin. In a later paper, 
Mullan & .lohnsorJ ( 199-5ll found oscillations in X-ray data of 


dMe stars. These oscillations have periods in the range of sev¬ 
eral tens to a few hundreds of seconds. They were not asso¬ 
ciated with any flare, though they were interpreted as coronal 
loop oscill ations. Optical stellar flare o scillations have been ob¬ 
served by iMathioudakis et all (l2003l) . who found a period of 
220 s in the decay phase of a white-light flare on the RS CVn 
binary II Peg. 

In this work, we investigate the first X-ray observation of 
an oscillation during a stellar flare. Section|2ldescribes the tar¬ 
get, AT Mic, as well as the observation. In Sect. 0 we present 
the data analysis, determining the period and amplitude of the 
oscillation through wavelet analysis. In Sect. |3 we determine 
the magnetic field strength and the length of the coronal loop 
by assuming that the oscillation is due to a magneto-acoustic 
wave. As a validity check, the value of loop length was com¬ 
pared to the value determined from a radiative cooling model, 
as well as from pressure balance scaling laws. The latter also 
gives an independent estimate of the magnetic field strength. 
Finally, we discuss the results and give conclusions (Sect.|3- 


2. Target and Observation 

2.1. Target 

AT Mic (GJ 799A/B) is an M-type binary dwarf, with both 
stars of the same spectral type (dM4.5e-i-dM4.5e). Both com¬ 
ponents of the binary flar e frequently. The radius of AT Mic 
given bv iLim e t alJlll98 7l) is 2.6 ■ 10*° cm . using a stellar dis¬ 
tance of 8.14 PC dOliese & .Tahreisslll99lll . Correcting for the 
newer value for the di stance from HIPPARCOS (10.2 + 0.5 pc, 
IPerrvman et aljfl99^ . and using the fact that the stellar radius 
is proportional to the stellar distance for a given luminosity and 
spectral class, we obtain a stellar radius of = 3.3 ■ 1 0^° cm = 
0.47 rQ. The mass of AT Mic given by iLim et all lll987l) is 
m* = 0.4 niQ. 

2.2. Observation 

For o ur analysis we used the XMM-Newton (l.lansen et alJ 
l200lh observati ons of AT M i c on 16 October 2000 during 
revolution 156. iRaassen et al.] il2003li have analysed this data 
spectroscopically, obtaining elemental abundances, tempera¬ 
tures, densities and emission measures, while a comparative 
flare analysis between X-ray and si multaneously observed ul - 
traviolet emissions can be found in iMitra-Kraev et ^ ( l2005l) . 
Here, we solely used the 0.2-12 keV X-ray 

S ’" iropean Photon Imaging Camera (EPIC-p 


data from the 
i.Lstriider et al.l 
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Fig. 1. The AT Mic 0.2-12 keV X-ray 10-s bin light curve 
(shown in black). Overplotted in white is the same light curve 
smoothed with a sliding time window of 200 s. The time is in 
ks, starting from the beginning of the observation (2000-10-16 
00:42:00). The vertical dotted lines mark the flare start, end of 
rise phase, end of extended top phase and end of the flare. 


The observation started at 00:42:00 and lasted for 25.1 ks 
(~7 h). Figure [2 shows the AT Mic light curve. A large flare, 
starting ~15 ks into the observation, increases the count-rate 
from flare onset to flare peak by a factor of 1.7, and lasts for 
1 h 25 min. It shows a steep rise (rise time = 1300 s) and 
decay (decay time - 1700 s). There is an extended peak to 
this flare, which shows cl ear oscillatory behavi our. Applying a 
multi-temperature model. iRaassen et al.lJ2003l) obtain a best fit 
with a mean flare temperature T = 24 + 4 MK and a quiescent 
temperature Tg - 13 + 1 MK, and from the O vii line ratio a flare 
and quiescent electron density of n = 4;|;3 10'° cm^^ and = 
(1.9+1.5) 10'° cm“^, respectively. The total flare and quiescent 
emission measures are EM = (19.5+0.8) 10^* cm~^ andEM^ = 
(12.2 + 0.5) 10^' cm^^. The total emitted energy in the 0-2- 
12 ke V band of this flare is ~ 6 ■ 10^^ erg dMitra-Kraev et al.l 
l200.5h . 

3. Analysis 

3.1. Data preparation 

Eor the extraction of the EPIC-pn light curve we used the 
XMM-Newton Science Analysis System (SAS) version 6.0. 
Light curves with a cadence of 1 s were extracted for con¬ 
centric regions around the source (20:41:51.156, -32:26:11.02, 
ICRS 2000 coordinates), and a background, off-set from the 
source on the same detector chip (20:41:48.428,-32:29:00.66); 
both regions used a radius of 41.325 arc sec. The selected en¬ 
ergy range is 0.2-12 keV (T = 1-62 A). The obtained errors 
for the count rates follow Poisson statistics for all data points. 
The background light curve is subtracted from the source light 
curve. Data missing due to telemetry absences are identified by 
obtaining the 1-s light curve of the entire detector chip; they 
manifest as time intervals with zero count rate. At such times, 
the background-subtracted raw light curve is interpolated. This 
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Fig. 2. The upper panel shows the 10-s bin light curve. The 
lower panel displays the high-frequency wavelet coefficients of 
the continuous wavelet transform of the above fight curve, us¬ 
ing a Morlet wavelet (see Sect. 1.3.2> . The contours give the 68%, 
95% and 99.9% significance levels. Tfie dasfied lines represent 
the cone of influence. The vertical dotted lines are the same 
as in Fig. ^ The arrow points to the flare oscillation in the 
wavelet domain and the horizontal dotted lines mark the di¬ 
vision between the high (noise), medium (oscillation) and low 
frequency ranges. 

fight curve is then re-binned to 10-s bins and used for further 
analysis. 

3.2. Continuous wavelet transform and frequency 
band decomposition 

We apply a continuous wavelet transform (CWT) to the 
10-s fight curve, follow ing the approach laid out by 
iTorrence & Compaill998l) . For our analysis, we use the Morlet 
wavelet Figure |2l displays the CWT. 

The vertical dotted fines denote the start (14.3 ks), the end 
of the rise phase (15.6 ks), the beginning of the decay phase 
(17.9 ks) and the end (19.6 ks) of the flare. The contours give 
significance levels of 68%, 95% and 99.9%, respectively. At 
flare rise, tfie power is enfianced at all timescales in tfie CWT. 
Tfie bold arrow points to the flare oscillation, which causes a 
local maximum in the wavelet coefficient plane. It has a pe¬ 
riod between 500 and 1200 s and occurs during flare maxi¬ 
mum, with a significance level of > 99.9%. Note tfiat tfiere 
is anotfier local maximum of 68% significance at a period of 
~1000 s occurring during a weaker flare at around f = 10 ks, 
wfiicfi suggests tfiat tfiis earlier flare migfit oscillate as well. A 
third, minor flare at f = 22 ks, on the other hand, shows no os¬ 
cillation. Here, we investigate only the major oscillation during 
the largest flare. 

We split the CWT into three frequency bands with periods 
of F > 1200 s, 500 s < P < 1200 s, and 10 s < F < 500 s, so 
that the middle band encompasses only the oscillation range. 
The fight curve is reconstructed for the three bands separately. 
Figure 0 shows these decomposed light curves (upper three 



Time [ks] 

Fig. 3. The reconstructed light curve divided into three fre¬ 
quency bands, low (F > 1200 s), medium (500 s < F < 1200 s, 
including the oscillation) and high (10 s < F < 500 s, mainly 
noise); together they add up to the original data (lowest panel). 
The shaded areas show the standard errors. The vertical dotted 
lines are the same as in Fig.[3 A coherent oscillation is easily 
identified during flare peak. 
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Fig. 4. Tfie flare part of tfie F = 500... 1200 s reconstructed ligfit 
curve (Fig. 0 second panel). Tfie vertical dotted lines bound 
the flare top. The solid fine sfiows tfie reconstructed fight curve 
from the data, whereas the dashed fine is a damped sine curve 
(bounded by the exponential envelope, dashed-dotted line) with 
an oscillation period of F = 750 s and an exponential decay 
time of T = 2000 s. The relative peak-to-peak amplitude is ini¬ 
tially 15%. 

panels), and their sum, the reconstructed fight curve (lowest 
panel). As the decomposition is a linear transform, it is easy 
to derive the standard errors for the decomposed light curves 
(see AppendixIXlfor the error treatment in CWT analysis), they 
are represented by the shaded areas around each fight curve. 
Figure 0 demonstrates how the oscillation has been separated 
from the flare profile, allowing the oscillation to be analysed 
separately. 
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Figure0]shows the isolated oscillation (solid line) to which 
we fit a damped sine curve (dashed line). The fitted curve has 
an oscillation with a period of P = 750 s, an exponential decay 
(damping) time of r = 2000 s, and an amplitude of 1 count s ', 
which gives, with an average count rate of 13 counts s“* during 
flare peak (from the low-frequency light curve), a peak-to-peak 
amplitude of A/// = 15%. 

A continuous wavelet transform is able to provide us with 
all our aims simultaneously: the wave is located in time and 
frequency (=l/period), it is naturally isolated from the underly¬ 
ing flare profile and thus the amplitude and damping are easily 
available from the fit parameters for a damped sine wave. The 
low signal-to-noise ratio is dealt with by disentangling the high 
frequencies from the lower oscillating frequency. In particular, 
this provides us with a better estimate for the amplitude. While 
the period derived from the smoothed light curve is broadly the 
same as the period obtained from wavelet analysis, the ampli¬ 
tude A/// — (7mc/.v 7mm)[(7mfl.Y "b 7mm)/2] — 10%, with Ifuax 

(Imin) the intensity at the first peak (dip), is considerably under¬ 
estimated. 

4. Results 

Our analysis clearly shows that there is an oscillation during the 
X-ray flare on AT Mic. The oscillation, which has a confidence 
level > 99.9%, starts when the flare reaches its maximum and 
continues during the extended flare top. A damped sine curve 
with an oscillation period of P = 750 s and an exponential 
decay time of t = 2000 s fits the data. The initial relative am¬ 
plitude is A/// = 0.15. 

In the following, we investigate the different types of 
magneto-acoustic loop oscillations. Comparing their various 
estimates for the flare loop length with two other, independent 
methods for the loop length, namely from radiative cooling 
times and from pressure balance, we find the most likely os¬ 
cillation mode. We also derive the magnetic field strength. 

4.1. Magneto-acoustic waves 

IZaitsev & Stenano^ ( Il989l) introduced a model, where a hard 
X-ray oscillation is excited by a centrifugal force, which the 
evaporating plasma exerts on the flare loop, stretching the mag¬ 
netic field lines upward and triggering Alfven oscillations. In 
this scenario, the relative amplitude of the oscillation is deter¬ 
mined by the additional energy from filling the magnetic flux 
tube with hot plasma 

AI AnnksT 

T “ p2 ■ ^ ^ 

With A/// = 0.15, n = (4 + 3) ■ 10'° cm^^^ P = (24 + 4) MK 
and ks the Boltzmann constant, the magnetic field then is 

(105+ 50) G _ 

[Roberts et alJ (Il984l) showed that for a straight cylindrical 
geometry several types of magneto-acoustic wave modes in a 
magnetic loop are feasible, namely the slow (acoustic), the fast 
kink and the fast sausage modes. For a standing oscillation in a 
loop, the loop length L is given by 

L = JcP/2, (2) 


with P the oscillation period, j the mode number and c the 
appropriate wave speed. For slow modes, c is the tube speed c,, 
with 



( 3 ) 


whereas for the fast kink modes, c is the kink speed Ck, with 


ncl + riecl^ 

n + He 


( 4 ) 


with the sound speed Cs - yfypTp and the Alfven speed 
CA = y - 5/3 is the adiabatic constant for a 

monatomic gas/plasma, p the gas pressure and p the plasma 
density. For the fast sausage mode, P = Inajck, where a is the 
tube radius. We set j - \ for the fundamental mode, where the 
loop apex has maxi mum displacement, in agre ement with the 
centrifugal model of IZaitsev & SteDanovllll989l) . 

The sound speed is Cs - ^jlyksTImp - 1.66T0"'Vp = 
(8.1+0.7) 10^ cm s“*. Using the centrifugal model and in¬ 
serting Eq. (□ into the expression for the Alfven speed, we 
obtain CA = (kBTy^^[mp(AI/I)]-^/^ = (1.1+0.1) 10^ cm s^'. 
We immediately see that, assuming the centrifugal model, the 
ratio of Alfven to sound speed is determined by the relative 
oscillation amplitude and is ca/cs = [ 27 (A ///)]“*^2 ~ 1.41. 
Using this last relation, the tube speed then is c, 0.82cs = 
(6.6 + 0.6) 10^ cm s *. To obtain the kink speed, we also need 
to know the Alfven speed outside the loop, cac- Assuming the 
same value for the magnetic field strength inside and outside 
the loop, CAe — ca ^rijne - (1.7 + 1.4) 10* cm s“*. Note the 
large error because of the large density error. The kink speed 
then is ct - ca yllnlin + «<,) = (1.3 + 1.1) 10* cms“*. Thus, 
we have (ct - 0.82cj) < Cs < (ca - 1.41Ci) < (ck = 1.65cj) < 
(cAe - 2.05Ci). 

Assuming a standing slow-mode oscillation, the loop 
length is Lsm - Pctl2 - (2.5 + 0.2) 10*" cm. For a standing 
fast kink-mode oscillation, a loop length of Lfk = Pckl2 - 
(5+4) 10*" cm is derived. And assuming a fast sausage 
mode oscillation, a - Pckli^n) - (1.6 + 1.4) 10*" cm and 
Lfs = (8...16) 10*" cm, assuming an aspect ratio of L/a = 5...10. 


4.2. Loop length from radiative cooling times 

The loop length can also (and independently of any oscillation) 
be estimated from rising and cooling times obtained from the 
temporal shape of the flare, applyi ng a flare heating/cooling 
m odel (see, e.g., lc argill et al.lll995l) . We follow the approach 
bv iHawlev et al.Nl995h who investigated a flare on AD Leo 
observed in the extreme ultraviolet. The shape of this flare is 
very similar to our flare on AT Mic, but roughly 10 times larger, 
and shows a flat top, too. 

The flare loop energy equation for the spatial average is 
given by 

\p-Q-R, ( 5 ) 

with Q the volumetric flare heating rate, R the optically thin 
cooling rate and p the time rate change of the loop pressure. 
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During the rise phase, strong evaporative heating is dominant 
(2 » R), while the decay phase is dominated by radiative cool¬ 
ing and strong condensation (R » Q). At the loop top, there is 
an equilibrium (R - Q). The loop length can be derived as 


L = 


1500 




4/7 


■4/7 . t-3/V .7-1/2 


(6) 


where Tr is the rise time, the flare decay time (indicated in 
Fig. |2]with the vertical dotted lines), T the apex flare temper¬ 
ature and = Cdicmax, with Cmax the peak count rate and Cd 
the count rate at the end of the flare. Inserting these values, the 
loop length is found to be L « 2.5 ■ 10*° cm. 


Table 1. Comparison of values derived for the indicated models 


Model 

Loop length 

Magnetic field 

Slow mode" 

(2.5 ± 0.2) 10'° cm 


Fast kink mode" 

(5 ± 4) 10'° cm 


Fast sausage mode" 
Centrifugal force model* 

(8...16) 10'° cm 

(105 ± 50) G 

Pressure balance^^ 

(2.8 ± 1.7) 10'° cm 

(75 ± 40) G 

Radiative cooling* 

2.5 • 10'° cm 



a 

Roberts et alJl984ll 

b 

Zaitsev & Stenanovl 19891) 

c 

Shibafa & Yokovamal2002h 

d 

Hawlev et alJl99a) 


4.3. Pressure balance 


To maintain stable flare loops, the gas pressure of the evapo¬ 
rated plasma must be smaller than the magnetic pressure 


2nkT <—. (7) 

Stt ^ 


Knowing the flare density and temperature, we get a lower limit 
for the magnetic field strength B > 80 + 60 G. Again, there is 
a large error because of the large uncertainty for the density. 
[^ib^& movain^ il^O^ assume pressure balance and give 
equations for B and L: 


B = 50 


L = 10** 


EM 


10481 


n-3 


EM 


J048| 


) ( 


tie 


\3/10 


/ T 


10° cm^ 

n. 


10 ° 


cm“ 


r[ 


17/10 
J x-8/5 

107 K 


\ 

) G (8) 

^-8/5 

j cm. (9) 


Using these relations, we obtain a magnetic field strength of 
B = 75 + 40 G and a loop length of L = (2.8 + 1.7) 10*° cm. 


Combining the above results, we find that both pressure 
balance considerations as well as the centrifugal oscillation 
model are consistent. The derived loop lengths from pressure 
balance and radiative cooling times are consistent with each 
other and in agreement with the loop length derived assuming 
a slow-mode oscillation or a fast kink mode. The fast sausage 
mode is much less likely to be the cause of this oscillation, 
as its derived loop length is inconsistent with the other inde¬ 
pendent methods. We further discuss the plausibility of the 
different modes in the following section. 


5. Discussion and conciusions 

We have investigated an X-ray oscillation during a flare on AT 
Mic, observed by XMM-Newton. In order to locate the oscilla¬ 
tion in time and frequency, we applied a continuous wavelet 
transform to the time series, and found that the oscillation, 
starting at flare peak and continuing during the flare’s flat-top 
phase, with a period around P = 750 s, has a confidence level 
of > 99.9% (Fig. I2J. We then reconstructed the light curve 
from the CWT for three different frequency bands (Fig. 0. 
In particular, the light curve passed through a low-pass filter 
(P > 1200 s, first panel), displays the overall shape of the 
light curve and picks out the long-lasting flares, whereas the 


light curve reconstructed from wavelet coefficients in the range 
around the oscillation period (500 s < P < 1200 s, second 
panel), isolates the oscillation. The coherent oscillation during 
the large flare can be fitted with a damped sine curve with a pe¬ 
riod of P = 750 s and an exponential decay time of r = 2000 s 
(Fig.0. The initial, relative peak-to-peak amplitude of the os¬ 
cillation is found to be AI/I = 15%. Decomposing a light curve 
in such a way is a powerful tool for isolating an oscillation from 
a flare, as well as disentangling the flare and the oscillation 
from the high-frequency noise. Because the effect of noise can 
thus be greatly reduced, this method provides us with a more 
accurate value for the oscillation amplitude than a value poten¬ 
tially obtained from a binned or smoothed light curve, which is 
always underestimated, especially if the bins are large. 

Interpreting this oscillation as a standing magneto-acoustic 
wave in the flare loop, we infer that it is a longitudinal slow¬ 
mode wave, oscillating at the fundamental frequency. This 
mode has an anti-node (maximum disturbance) at the loop apex 
and is capable of causing th e largest global flux variations. 
Using the relations derived bv [Roberts et al.l lll984l) . we find a 
flare loop length of L = (2.5 + 0.2) 10*° cm. This value is con¬ 
sistent with estim ating the loop leng^ th from radiative cooling 
times of the flare ilFfawlev e t al.ll 99 ^ as well as from pressure 
balance considerations dShihata & YokovamalboO^ . Table ^ 
provides a comparison of these results. The derived loop length 
of a fast kink wave, which has a large error due to a large error 
in density, is also consistent with the two independent meth¬ 
ods of radiative cooling times and pres sure balance. This mode , 
however, is basically incompressible dNakariakov et al.ll2004l) . 
and it would be hard to imagine how it could cause intensity 
perturbations in a spatially unresolved light curve. We discount 
fast sausage-mode waves, since they would require either much 
longer loops (several times the stellar radius) for the observed 
period, which is inconsistent with the loop length derived from 
the radiative cooling as well as from pressure balance mod¬ 
els, or a loop aspect ratio (L/a) close to 1, in which case the 
cylindrical geometry assumption breaks down. In summary, we 
conclude that the observed oscillation is most likely a standing 
longitudinal slow-mode wave. 

To estimate the ma gnetic field strength, we app ly the cen¬ 
trifugal force model of IZaitsev & StepanovI lll989t) . where the 
magnetic loop is stretched by a centrifugal force which is 
caused by the upwardly evaporating plasma, stretching the loop 
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beyond equilibrium and thus exciting MHD waves. The flux 
tube starts to oscillate up and down (i.e. in the plane of the 
loop). Consequently, the particle density inside the loop also 
oscillates, causing an oscillation in the thermal radiation (soft 
X-rays and extreme ultra-violet). Non-thermal radiation, origi¬ 
nating from gyro-synchrotron (radio) and bremsstrahlung (hard 
X-rays and optical), caused by fast electrons trapped within 
the loop and moving back and forth from one end of the loop 
to the other, also oscillates: the speed of the trapped elec¬ 
trons, and hence the non-thermal radiation, is modulated by 
the plasma density. This picture is consistent with a slow-mode 
(acoustic) wave. The magnetic field strength thus obtained is 
B - (105 + 50) G. The large error results from the large un¬ 
certainty in the particle density. This value for the magnetic 
field strength is consistent w ith pressure balance considerations 
dShibata & Yokovamal2002l) . where B = (75 + 40) G. 

The observed rapid damping of the oscillation (a damp¬ 
ing time comparable to the oscillation period) is in agree- 
ment with other magneto-acoustic wave ob servations (e.g.. 


ment witn otner magneto-acoustic wave ob servation 
Ofman & Wand 2002t IVerwichte et alJ 20041) . So far. 


models have been put forwar d to explain the fast dam ping. In 
the case of slow-mode waves, lOfman & Wand (l2002h numeri¬ 
cally find that thermal conduction is the dominant dissipation 
mechanism in loops with T > 6 MK. 

An alternative to the interpretation of a standing magneto¬ 
acoustic wave causing the flare oscillation is that repeated and 
rapid flaring is occurring. This might also explain the flat-top 
character of the time profile of this flare, which differs from 
the usual shape of solar flares, which shows a rapid r ise from 
begin ning to peak, followed by a slow decay (see, e.g., ISvestkal 
19891) . However, repeated flaring should be random, rather than 


have a periodic and exponentially-decaying nature, and would 
have difficulties explaining the observation of what appears to 
be a damped oscillation. 

The derived length of the flare loop is about the size of the 
stellar radius and is c omparab l e to l oop lengths of large so¬ 
lar flares. In narticular. lSvestkal ill 9941) investigated solar X-ray 
observations which showed periods close to 20 min. They as¬ 
sociated them with spatially resolved large-scale coronal loops 
with loop lengths of (2-3) 10'° cm, and interpreted t he os- 
cillations as slow-mode MHD waves. iTerekhov et al ] (I2OO2I) 
observed an X-ray oscillation with P - 143 s during a so¬ 
lar flare, and derived a loop length of (1-3) 10*° cm. The 
centrifu gal force model, originall y derived for non-thermal 
X-rays jZaitsev & Stepan^ 1 19891) , has also bee n applied to 
milli metre-wave emission dstenanov et alJ[T992t) . thermal X- 
rays ijTgr^hov^t^ ll2Q02l). and, in the stellar case, to optical 
emission dMiillan et al.ll9^[Mathioudakis et al.i2003l). 

This is the first time that an oscillation has been observed 
in X-rays in a stellar flare, and has been used to derive a flare 
loop length and a local magnetic field strength. Comparable 
loop dimensions and magnetic field strengths indicate the sim¬ 
ilar natures of AT Mic’s corona and that of our Sun. Because 
of the smaller radius of these M-type stars (r* » 0.5 tq), the 
loop length is in fact of the order of the stellar radius, whereas 
the magnetic field strength is around the upper limit for solar 
coronal values. This is consistent with this class of low mass 
stars (dMe-type) being very X-ray active. 
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Appendix A: CWT error analysis 

The continuous wavelet transform of the signal / is given by 
convolving the signal with a wavelet function 'P 

(A.l) 


or in components: 

(A.2) 

m 


with Wjn = W{tn,Sj), 'Py,™ = 

fm - f{tm), the * denoting the complex conjugate and the scale 
Sj = So In our analysis, we set sq = 1 and Sj = 0.25 for 
the Morlet wavelet. Significance levels can be derived as 



(A.3) 

(A.4) 

(A.5) 


for cr^ the variance of /,„, and the -distribution of Wjn- 
The signal can be reconstructed by using 


y-l P.e{W jn] 
/« = C ■ > -— 


(A.6) 


with c the norma lisation constant (see, e.g., 
iTorrence & Compel Il998l) . If a filter F is applied to the 
wavelet coefficients, Wjn - Wjn ■ Fjn, the filter coefficients 
being 0 or 1, the filtered signal is then given by 


'Re{Wjn} 


(A.7) 


Inserting Eq. ) IA.2t into Eq. and rearranging the sums, we 

obtain 


fn — ^ ] ^mn fm 


(A.8) 


with 


Amn — ■ 'y 




jnm i ‘ Fjn 




(A.9) 


Thus, the signal can be filtered and reconstructed using a linear 
transform, which is determined by the wavelet function and the 
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filter. The variance of the reconstructed light curve is then given 
by 

m 

The signal can be decomposed without loss, if the sum of the 
different filters used is 1 for all filter components. 
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